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Abstract. Spectral tetris is a flexible and elementary method to construct unit norm 
frames with a given frame operator, having all of its eigenvalues greater than or equal to 
two. One important application of spectral tetris is the construction of fusion frames. We 
first show how the assumption on the spectrum of the frame operator can be dropped and 
, ^ . extend the spectral tetris algorithm to construct unit norm frames with any given spectrum 

of the frame operator. We then provide a sufficient condition for using this generalization 
of spectral tetris to construct fusion frames with prescribed spectrum for the fusion frame 
operator and with prescribed dimensions for the subspaces. This condition is shown to be 
, necessary in the tight case of redundancy greater than two. 

1. Introduction 

A fusion frame is a sequence of subspaces of a Hilbert space along with a sequence of 
weights, so that the sequence of weighted orthogonal projections onto these subspaces sums 
to an invertible operator on the space. Fusion frames, introduced in [7] and refined in [9], have 



O 



become a subject of interest due to their applicability to problems in distributed processing, 
sensor networks and a host of other directions. Fusion frames provide resilience to noise 
and erasures due to, for instance, sensor failures or buffer overflows [TJ El ETUI EE], as well as 
00 ■ robustness to subspace perturbations [9] , which can happen because of imprecise knowledge 
of sensor network topology. For fusion frame applications, we generally need extra structure 
on the fusion frame, such as prescribing the fusion frame operator or the dimensions of the 
subspaces, or both. 

In this paper we address the question of how to efficiently construct fusion frames with 
prescribed dimensions of the subspaces and prescribed eigenvalues of the fusion frame oper- 
ator. Our main tool is the spectral tetris construction for unit norm frames, which we will 
review in section This construction is limited to the case of frames having a frame operator 
with spectrum in [2, oo) and therefore we will first, in section [31 extend spectral tetris to the 
most general possible case in terms of the prescribed spectrum of the frame operator of the 
unit norm frame. Precisely, we give a version of spectral tetris, capable of constructing an 
M-element unit norm frame in with prescribed eigenvalues (A n )^ =1 C (0, oo) satisfying 
only the necessary trace condition J2n=i ^« = M. Before tackling this case in section 13.21 we 
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first consider the case of tight frames in section 13.14 i.e. we extend the existing construction 
to unit norm tight frames of redundancy less than two. 

In section HI we give a necessary condition on a prescribed sequence of eigenvalues and 
dimensions, under which we can use the generalized spectral tetris method to construct a 
fusion frame having those eigenvalues for its fusion frame operator and having those dimen- 
sions for its subspaces. We further show that this condition is also necessary in the case of 
unit norm tight frames of redundancy of at least two. 



2. Preliminaries and notation 

2.1. Fusion frames. Let (e n )^ =l be the standard unit vector basis of K . The synthesis 
operator of a finite sequence (fm)m=i — ^ N * s ^ ' ^ ~~ ^ ^ N & ven by 

M 

Fg = ^2(g,e m )f m , 

771=1 

i.e. F is the N x M matrix whose m-th column is f m . The sequence (f m )m=i ls a / rame if 
its frame operator S = FF* satisfies AI < S < BI for some positive constants A, B, where 
/ is the identity on <C N . In particular, the spectrum of S is positive real. The sequence is a 
tight frame if A = B, i.e. if 

M 

Af = / ]{fi fm)fm 

771=1 

for all / G C^, or equivalently if 

M ( A I 

S^/f \77 \-J ' ra = n ' 

^— ' 0, n . 

771=1 ^ ' 

In the case of a tight frame, the constant A equals M/N and is also called the tight frame 
bound or the redundancy of the frame. The synthesis matrix of the frame (f m )^ =1 is called 
s-sparse, if it has s nonzero entries. A unit norm tight frame is a tight frame {fm)m=i 
for which ||/ m || = 1 for all m = 1,...,M. Unit norm tight frames provide Parseval-like 
decompositions in terms of nonorthogonal vectors of unit norm. If {fm)m=i * s a um ^ norm 
frame, the operators / \-> (f, f m )f m arising in the frame operator 

M 

Sf = ^ fm)fm 

771=1 

are rank-one orthogonal projections. Fusion frame theory is the study of sums of projec- 
tions with weights and of arbitrary rank. In particular, a sequence (Wfe, ffc)fc=i °f subspaces 
{Wk)k=\ °f C N and weights (v k )f =1 is a fusion frame if the sequence {Pk)k=\ °f orthogonal 
projections onto those subspaces satisfies 



K 

AI < J2 v l p k < BI 

k=l 
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for some positive constants A, B. It is a tight fusion frame if A = B. Using the Horn- 
Klyachko compatibility inequalities, [TT] gives a characterization of the sequences of weights 
and dimensions of the subspaces for which tight fusion frames exist. This result however can 
hardly be used in practice, since it involves the computation of the Littlewood-Richardson 
coefficients of certain associated partitions and moreover does not give a construction of the 
existing fusion frames. In this paper we restrict ourselves to the case where all weights are 
equal to one and denote the fusion frame by (Wk)f =v In this case, the fusion frame operator 

is S = X^fcLi Pk- If Dk is the dimension of the range of Pk and (fk,d)d=i is an orthonormal 
basis of the range of Pk then 

K K D fe 

sf = y~]Pkf = y^ J 'y^ J {f,fk < d)fk,d 

k=l k=l d=l 

for all / G C^. This shows that every fusion frame arises from a classical frame that 
satisfies additional orthogonality requirements. To be precise, we say that a sequence 
{fk,d)k=id=i — ^ N generates a fusion frame (Wk)% =1 with dimPVfc = D k for k — 1 . . . , K 
L =1 d= \ is a frame for C N and (fk,d)d=i i s orthonormal for k — 1 . . . , K. 
We say that a (fusion) frame has a certain spectrum or certain eigenvalues, if its (fusion) 
frame operator has this spectrum, respectively these eigenvalues. Note that we are always 
counting multiplicities of the eigenvalues. For any frame, the sum of its eigenvalues equals 
the sum of the lengths of its vectors, i.e. the number of vectors in case we are dealing with 
unit norm frames. 

Building on the theory of majorization and the Schur-Horn Theorem, [2] shows how to 
explicitly construct every possible frame whose frame operator has a given arbitrary spectrum 
and whose vectors are of given arbitrary lengths. In this paper we use a much easier algorithm 
for the construction of unit norm frames with given spectrum. This algorithm is called 
spectral tetris and we will discuss it in detail below. We focus our attention on this algorithm, 
since it is tailor made for the construction of fusion frames, as it constructs frames of vectors, 
a lot of which are orthogonal to each other due to disjointness of their supports. 

2.2. Spectral Tetris. The term spectral tetris refers to the first systematic method for 
constructing unit norm tight frames. This construction was introduced in [6] to generate unit 
norm tight frames in M. N for any dimension N and any number of frame vectors M, provided 
that M > 2N. Choosing all weights to equal 1, [6] provides a complete characterization of 
triples (N, K, d) for which tight fusion frames of K subspaces of equal dimension d exist in 
R and gives an elegant algorithm to produce such tight fusion frames for most of those 
triples. 

A straightforward extension to the construction of unit norm frames having a desired frame 
operator with eigenvalues (A n )^ =1 C [2, oo) satisfying J2n=i^ n = M was introduced in [3]. 
For convenience we review this construction in table [Hand will refer to it as the spectral tetris 
construction (STC). STC constructs synthesis matrices with unit norm columns, whose rows 
are pairwise orthogonal and square sum to the desired eigenvalues. This way, the frame 
operator is a diagonal matrix, having precisely the desired eigenvalues on its diagonal. We 
strongly recommend a glance at [6] or [3] for instructive examples on how the algorithm 
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constructs the desired synthesis matrices by using 2x2 building blocks of the form 



A construction for equi-dimensional fusion frames having eigenvalues as above for their fusion 
frame operators is given in j3] . The sufficient condition for this construction to work, is that 
the dimension d of the subspaces satisfies ^2 n=1 ^« = dK where K is the number of subspaces 
and A n < K — 3 for n — 1, . . . , N. 



STC: Spectral Tetris Construction 

Parameters: 

• Dimension N G N. 

• Number of frame elements M e N. 

• Eigenvalues (A n )^ =1 C [2, oo) such that ^2 n=1 A n = M. 
Algorithm: 

1) Set k = 1. 

2) For j = 1, . . . , N do 

3) Repeat 

4 If A,- < 1 then 



(1) 



\/l — x —VI — x 



6 



5 




7 

8 
9 

io : 

11 : 
12 : 
13 
14 

15 : 
16 : 



k = k + 2. 



Aj+i — Aj + i — (2 
A, = 0. 




else 



fk = ej- 
k = k + l. 

= ~ !• 



end 
until A, = 



0. 



end. 



Output: 



• Unit norm frame (fk)kLi Q K 



TABLE 1. The STC algorithm for constructing a unit norm frame with pre- 
scribed spectrum in [2, oo). 



In this paper, we extend the spectral tetris construction to the case of arbitrary eigenvalues 
for the frame operator and use those frames for the construction of fusion frames with given 
fusion frame operators and subspaces of not necessarily equal dimensions. 
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Definition 2.1. A frame constructed via the spectral tetris construction STC is called a 
spectral tetris frame. A fusion frame {Wk)f =1 is called a spectral tetris fusion frame if there 

is a partition of a spectral tetris frame (fk,d)k=id=i suc h that (fk,d)d=i ^ s an orthonormal 
basis for for every k = 1, . . . , K . 

Aside from the fact that spectral tetris frames are easy to construct, their major advantage 
for applications is the sparsity of their synthesis matrices. This sparsity is dependent on the 
ordering of the given sequence of eigenvalues for which STC is performed. Note that the 
original form of the algorithm in [3] assumes the sequence of eigenvalues to be in decreasing 
order. However, this assumption was made only for classification reasons, and it is easily 
seen that it can be dropped. The sparsest synthesis matrices are achieved if the sequence 
of eigenvalues (A n )^ =1 is ordered blockwise, i.e. if for any permutation ir of {1, . . . , N} the 
set of partial sums {X^=i '■ s = !>•••> N} contains at least as many integers as the set 
{^2 S j=i : s = • • • ' ^ nas been shown in [5], that spectral tetris frames are optimally 
sparse in the sense that given M > 2N and a sequence of eigenvalues (A n )^ =1 C [2, oo), the 
synthesis matrix of the spectral tetris frame having these parameters is sparsest in the class 
of all synthesis matrices of unit norm frames that have these parameters, provided STC is 
run for the sequence (A n )^ =1 rearranged to be ordered blockwise. Note that for tight frames 
all eigenvalues are equal, and as such, questions of rearranging the order of the eigenvalues 
do not arise. 

3. Spectral tetris for arbitrary prescribed spectra 

The spectral tetris construction STC discussed in section [2] is capable of constructing unit 
norm frames, all of whose eigenvalues are at least 2. It constructs the synthesis matrices of 
such frames by successively filling an appropriate sized matrix with ones and 2x2 blocks 
of the form ([I]). If we drop the condition on the spectrum to be contained in [2,oo), we 
can in general no longer use just 2x2 blocks for a spectral-tetris-like construction. We will 
instead work with larger building blocks, constructed from appropriate sized discrete Fourier 
transform matrices. 

Definition 3.1. Given L e N, let u = exp(^p) be a primitive L-th root of unity. The 
(non-normalized) discrete Fourier transform (DFT) matrix in C LxL is defined by 

Let us point out that here we do not normalize Fl by the factor 1 / a/Z, thus every entry 
of Fl is of modulus one. The importance of DFT matrices to our construction stems from 
the fact that they have orthogonal rows and columns, and that all entries have the same 
modulus. In the course of the construction we will have to alter the row norms of the DFT 
matrices by multiplying rows with appropriate constants. While this will destroy the pairwise 
orthogonality of the columns, it will preserve the pairwise orthogonality of the rows, which 
is the crucial feature for our construction to work. 

In all the constructions that follow, whenever a 2 x 2 DFT matrix is used, one might as 
well use a 2 x 2 matrix of the form ([I]). It is however not obvious how one could work with 
real matrices of size larger than 2x2 instead of complex DFT matrices. 
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It is perhaps most instructive to first look at an example of the construction we are going 
to introduce in this section. 



Example 3.2. We construct a 5 -element unit norm tight frame in C 4 . In the following we 
will use the notation ul = exp(^p). 

We can start filling the desired 4x5 synthesis matrix with an altered 2x2 DFT matrix 
in the upper left corner. The alteration we make, is to multiply the entries of the first row 
by a/5/8 in order to make the first row have the desired norm \/5/4. To get normalized 
columns, we multiply the second row of the 2x2 DFT matrix by a/3/8. At this point we 
have constructed the first row and the first two columns of the desired synthesis matrix: 







• 0J2 













Note that so far we have constructed a matrix whose first two rows are orthonormal, no 
matter how we keep filling the second row. The second row at this point has norm a/3/4, 
while we need to make it have norm a/574. We can not insert another alterated 2x2 DFT 
matrix in the same fashion as above, since we would have to multiply its first row by the 
factor a/2/8 in order to have the second row of the synthesis matrix square sum to 5/4, and 
its second row by the factor a/3/4 to get normalized columns. But then the second row of 
the altered DFT block would already square sum to 3/2 and thus exceed what we desire for 
the row norms of our synthesis matrix. We can however, alter a 3 x 3 DFT matrix in the 
above fashion. We multiply its first column by a/1/6 to add the missing weight to the second 
row of the desired synthesis matrix and make it square sum to 5/4. Moreover, we multiply 
the second and the third row of the 3x3 DFT matrix by a/5/12 to get 








• w 2 






19 \ 19 " U 3 \ 79 • ^, 



777 • 



The latter matrix is the synthesis matrix of the desired frame, since its columns are normal- 
ized and its rows are pairwise orthogonal and square sum to 5/4. 



We would like to extend definition 12.11 to include frames like the one in example I3.2[ i.e. 
in addition we call all frames that are constructed using altered DFT matrices as building 
blocks of their synthesis matrices spectral tetris frames, and fusion frames constructed by 
using the vectors of any such spectral tetris frames spectral tetris fusion frames. 

3.1. Tight low redundancy spectral tetris frames. We now give a version of spectral 
tetris for the construction of unit norm tight frames of redundancy smaller than 2, i.e. our 
goal is to construct an M-element unit norm tight frame for C N where N < M < 2N. In 
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most cases, the original construction STC does no longer work if M < 2N. Indeed, it is 
shown in [I], that in case N < M < 2N and M and N relatively prime, STC can construct 
a unit norm tight frame of M vectors in N dimensions, if and only if M — 2N — 1. Since, it 
is thus no longer possible to always work with 2x2 building blocks and ones, we will first 
address the question: What size alterated DFT-matrices should we use as building blocks of 
the to be constructed synthesis operator, in order to achieve the best possible sparsity? We 
want to develop an algorithm following the lines of example I3.2[ i.e. we want to compose the 
desired synthesis matrix out of square blocks whose first and last rows successively overlap. 
If we use K square matrices as building blocks of the N x M-synthesis matrix, then due to 
the overlapping, we have 

M = N+ (K- 1). 

We will first check the best possible sparsity we can get from this set-up, without considering 
at this moment the question, whether we can alter the DFT blocks in such a way, that the 
row norms of the synthesis matrix equal M/N and that the columns are normalized. 
Note that, if we do a spectral-tetris-like construction as in example 13.21 with K altered DFT 
matrices of sizes mi, . . . , trk, then the sparsity of the synthesis matrix we get is Ylk=i m fc- 

Lemma 3.3. Let m x < • • ■ < m^, and M be positive integers. The minimizer of Y^k=i m 1> 
subject to the constraints Ylk=i m k = and 1 < < K for k = 1, . . . , K , is given by 

mi — ■ ■ ■ — m r — L, and m r+1 = ■ ■ ■ = m K — L + 1, 
where L = [^\ and r = K(L + 1) - M. 
Proof. We first show that the minimizer above satisfies 

m K — m i + 1. (2) 
Indeed, assume we have achieved the minimum sum above and mx — mi > 2. Then 
(mi + l) 2 + (m^ - l) 2 = mj + m 2 K - 2(m k - m 1 - 1) 

< ml + m 2 K . 

Therefore, replacing mi and m^ by mi + 1 and m^ — 1, respectively, we would get a set of 
numbers satisfying the constraints of the lemma but having smaller square sum. Thus we 
have shown (J2]). 

Successive application of (J2J) implies, that there is an L and an index 1 < r < K, so that 
mi = • • ■ = m r = L, and m r+ i = • • • = mx = L + 1. 

Hence, 

K 

M = ^2m r = rL + {K - r){L + 1) = K(L + 1) - r. 

r=l 

Thus, r = K(L + 1) - M and since 1 < r < K, it follows that L = [f J. □ 

While STC used building blocks of sizes lxl (namely ones) and 2x2, lemma I3T31 suggests 
that we work with DFT matrices of the sizes L x L and (L + 1) x (L + 1) as building blocks 
to construct sparse unit norm tight frames of redundancy smaller than two. Here are the 
alterations we have to perform on those DFT matrices. 
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Definition 3.4. Let N < M < 2N and L = [ M ^ +1 \ ■ We denote a matrix by D L , or call 
it a Dl block, if it is derived from the DFT matrix Fl by multiplying the entries of the j-th 

row of Fl by J jjL, where c 2 = ■ • • = Cl-i = M and 1 < c±, Cl < M are integers, and if it 

has normalized columns. We call c\ and cl the first, respectively the last correction factor 
of Dl. In the same way, replacing L by L + 1, we define Dl+i, i.e. we multiply the j-th row 

°f F L+i ^ ^n(lW)- 

Note that the row norms of a D L or D L +i block equal M/N, except possibly for their 
first and last row. We need to collect some information about the first and last correction 
factors. 

Lemma 3.5. Let N < M < 2N and L = [ M M N+1 \ ■ 

(i) An integer ci is a first correction factor of a Dl block, if and only if L(N — M) + M < 
ci < M. 

(ii) IfL(N-M) + N>0, then any L(N-M) + N < c 1 < M is a first correction factor of 
aD L+1 block. IfL{N-M) + N < 0, then any L(N — M) + N < a < L(N-M) + M 
is a first correction factor of a D^ +1 block. 

Proof, (i) Let c\ be the first and cl the last correction factor of a block. We need to 
show 1 < cl < M. Since Dl has normalized columns, we have 

c L = NL-(L- 2)M - ci > NL - (L - 1)M = L(N — M) + M > 0, 

where the last inequality holds since L < . On the other hand, cl < M implies 

NL — (L — 2)M - ci < M, 

and thus L(N - M) + M < c x . 

(ii) Let Ci be the first and c^+i the last correction factor of the Dl+i block. The condition 
that the columns are normalized, is equivalent to 

CL+1 = N(L + 1) - (L-1)M - Cl . 

Thus cl+i < M is equivalent to 

N(L + 1) - (L - 1)M - ci < M, 

i.e. 

L{N -M) + N < ci. 

If L(N - M) + N > 0, then 

c L+1 = N(L + 1) - (L - 1)M - ci > N(L + 1) - (L - l)M - M = L(N - M) + N > 0. 
Now assume L(N - M) + N < and c x < L(N - M) + M. Then 

c L+1 = N(L + 1) - (L - l)M - ci > N(L + 1) - (L - 1)M - L(N - M) - M = N > 0. 

□ 

Definition 3.6. Let c±, cl 1 be the first, respectively last correction factor of a Dl 1 block and 
di,dL 2 be the first, respectively last correction factor of a Dl 2 block. If cl 1 + d\ = M then 
we call ci — di the stepsize of the D Ll block. 
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Lemma 3.7. Let N < M < 2N and L = [ M _ A ^ r+1 J • The step size of a Di block is 
L(N - M) + M, while the stepsize of a D L+1 block is L(N - M) + N . 

Proof. Consider first the case of a block with first correction factor c\ and last correction 
factor cl. By the normalization of the columns we have cl = LN — c\ — (L — 2)M . So the 
step size is 

ci - (M - c L ) = L(N -M) + M. 
In the case of a Dl+i block with correction factors C\ and Cl+i, we have by the column 
normalization cl+i = N(L + 1) — c\ — (L — 1)M, and thus for the step size 

ci - (M - c L ) = L(N -M) + N. 

□ 

In case [ M _^ +1 J = 1 we want to use the above mentioned result from [1], which charac- 
terizes the cases in which STC works, given that N < M < 2N. To this end, we make the 
following observation. 

Lemma 3.8. Let N < M < 2N be positive integers. Then M = 2N — 1, if and only if 
I — I = i 

Proof. Let 1 < k < N - 1 such that M = 2N - k. If k = 1, then 

M 

_M-N + l 

If k > 2, then 

M 

M-N+l 

□ 

We are now ready for the main result of this section. Given N < M < 2N, we propose 
the algorithm TDFTST (table [2]) to construct an M-element unit norm tight frame in C N . 
We analyze this algorithm in the proof of the following theorem. 

Theorem 3.9. Let M,N be relatively prime and N < M < 2N. Let K = M - N + 1, 
L = |_t^J and r = K(L + 1) — M. Then the synthesis matrix of the sparsest M-element 
spectral tetris unit norm tight frame in C N consists of r blocks and (K — r) blocks -Dl+i 
and can be constructed by DFTSTC. Its sparsity is rL 2 + (K — r)(L + l) 2 . 

Proof. If L — 1, the algorithm uses STC to construct the sparsest unit norm tight frame. 
This is possible since by lemma [378| L = 1 is equivalent to M — 2N — 1, which by [1] is the 
characterization of the cases STC can be used, given that M < 2N. 

Now suppose L > 2. In this case the algorithm builds the synthesis matrix of the desired 
frame as in example 13.21 by inserting and -Dl+i blocks. The last row of a block and the 
first row of the following block appear in the same row of the synthesis matrix and, while the 
last correction factor is chosen to ensure normalized columns, the following first correction 
factor is chosen to guarantee that the rows of the synthesis matrix square sum to M/N. As 
the algorithm progresses, we keep track of the development of the first correction factors in 
the variable x. We have to ensure that, whenever the algorithm is determined to insert a 





2iV-l 




N 



= 1. 



2N- k 
N-(k-l) 



> 2. 
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TDFTST: Tight DFT Spectral Tetris 


Parameters: 






• L^imeilSlOIl IV t 11 






• Number of frame elements 


l\l 


t wnere iv < ivi < ziv. 


- r 1 M 1 t jy- _ jf 

• Li — |_7?"J ' wiieie t\ — ivi — 


AT 
1 V 


i i 

T !■ 


Algorithm: 






1) If L = 1 






2) run STC with A n = § 


for 


n = l,...,N. 


3) else 






4) If L(N - M) + N < 






5) run subroutine 1. 






6) else 






7) run subroutine 2. 






8) end. 






9) end. 






Output: 






• Unit norm tight frame (f m 


\M 
)m= 


=1 C 



TABLE 2. The TDFTST algorithm for constructing a unit norm tight frame 
of redundancy less than two. 



certain block, the variable x that we computed in order to make the rows of the synthesis 
matrix square sum to M/N, lies in the range of integers that are possible as first correction 
factor for a block of the desired size, see lemma 13.51 The algorithm starts by letting the 
first correction factor x of the first block to be inserted into the synthesis matrix equal M. 
The difference between a first correction factor used in the course of the algorithm and the 
subsequent first correction factor that is being used, is the step size determined in lemma 
13.71 We have to ensure that eventually the sum of the step sizes equals M, i.e. the last 
correction factor of the final block inserted into the synthesis matrix equals M (in other 
words, the next first correction factor would be zero, but we have arrived at the point were 
the algorithm terminates). 

Let us begin by recording that the step sizes of r blocks Di and (K — r) blocks -Dl+i add 
up to M. Indeed, 

r(L(N -M)+M) + (K - r)(L(N - M) + N) 

= (KL - r )(N-M) + KN 

= (KL — K(L + 1) + M)(N — M) + KN (3) 

= M(N -M + K) 

= M. 
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Subroutine 1: 

1) m = 0,n=l,x = M 

2) Repeat 

3) If x > L(N -M) + M then 

4) u) = exp(^f) 

5) For r = m + I, . . . ,m + L do 

6) f - ,/~^e 4- * /-^- V" +L " 2 M0'-")(''-m-i)p. 

U J Jr — \J NL ^n ^ y NL l^j=n+l M e J 
j NL~M(L~2)~x (L-l)(r-m-l) e 

7) end. 

8) m = m + L 

9) n = n + L - 1 

10) x = x-L{N-M)-M 

11) else 

12) w = exp(|^) 

13) For r = m + l,...,m + L + l do 

14^1 f - I x r 4- / ^ V W+L ~ 1 LJ (.7-wXr-m-lU 

1 ^ J 7 " A/ JV(L+1) n V JVYL+1) Z_/j=n+l e J 



jV(L+l)-M(L-l)-:r , r( r _ m _!) 
AT(L+1) W 



15) end. 

16) m = m 4- L + 1 

17) n = n + L 

18) x = x-L{N-M) -JV 

19) end. 

20) until x = 0. 



Table 3. Subroutine 1 used in the TDFTST algorithm. 



To show the above mentioned properties, we now separately consider the cases of L(N — 
M) + N being positive and negative. 

We begin with the case L(N — M) + N < 0: In this case the algorithm first inserts a 
Dl block with first correction factor x = M. It then tracks the development of the first 
correction factors by subtracting the step size L(N — M) + M, respectively L(N — M) + N, 
whenever a Dl, respectively Dl+i block, has been inserted. The algorithm inserts a Dl 
block whenever possible, that is, whenever the first correction factor at a certain run of the 
repeat loop is at least L(N — M) + M. If and only if during a run of the repeat loop the 
first correction factor is less then L(N — M) + M, the algorithm inserts a Dl+i block and 
subtracts the step size L(N — M) + N from the first correction factor in order to get the first 
correction factor for the next run of the repeat loop. Note that, if the correction factor x 
falls below L(N — M) + M, it will still be positive and thus can serve as first correction factor 
for a D L +i block, since by lemma [33] it suffices that L(N - M) + N < x < L(N - M) + M 
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Subroutine 2: 

1) m = 0, n = 1, z = 

2) Repeat 

3) u = exp(^p) 

4) For r — m + 1, . . . , m + L do 

5) f - /~^ e + . [M. \^ n+L - 2 rjU-nKr-m-l) 
°) Jr — V NL^n -T \ NL l^j=n+l U e J 



6) end. 

7) m = m + L 

8) n = ra + L- l 

9) x = x- L(N -M)-M 

10) ^ = ^ + 1 

11) until z = K(L + l) - M 

12) Repeat 

13) u = exp(g) 

14) For r = m + l,...,m + L + l do 

151 f = I x e — t— a I M y >n+z '- 1 / < .(7-w)(r-m-i) e . 

1J / ^ r y Af(L+l) e ™ ^ y JV(i+l) 2^j=n+l W e J 
/ iV(L+l)-AJ(L-l)-a L(r-m-l) p 

+ Y W(lTT) U e n+L 

16) end. 

17) m = m + L + 1 

18) n = n + L 

19) x = x - L(iV - M) - iV 

20) z = z + 1 

21) until z = K. 

Table 4. Subroutine 1 used in the TDFTST algorithm. 



and since L(N - M) + N < 0. Since L(N - M) + N < 0, according to ©, the algorithm 
will thus terminate after inserting r blocks Dl and (JiT — r) blocks -D^+i in some order. 

Now consider the case L(N — M) + N > 0: In this case the algorithm fills the synthesis 
matrix by first using r blocks Dl followed by (K — r) blocks -Dl+i, in this order. By ([3]) we 
have 

M - r(L(N -M)+M) = (K- r)(L(N - M) + N) > 0, 

and thus 

M — (r — l)(L(N - M) + M)> L(N — M) + M, 

which implies that starting with correction factor x = M and inserting (r — 1) blocks Dl, 
i.e. subtracting (r — 1) step sizes L(N — M) + M from x, results in an integer greater than 
or equal to L(N — M) + M, which therefore can serve as the first correction factor of the 
last Dl block that is being inserted by the algorithm. After inserting this last Dl block, 
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(K — r) blocks Dl+i are being inserted. By ([3]), the final first correction factor prior to the 
termination of the algorithm is L(N — M) + N, i.e the final last correction factor is M. Thus 
the algorithm succeeds in constructing the desired synthesis matrix. □ 

Remark 3.10. If we drop the assumption on the number M of vectors and the dimension N 
of the space to be relatively prime, we can easily construct a unit norm tight frame by using 
gcd(M, N) copies of a matrix that TDFTST produces. To be precise, given positive integers 
N < M < 2N, let M' = M/gcd(M, N) and N' = N/gcd(M,N). Let F be the N' x M' 
synthesis matrix that TDFTST generates for M' vectors in N' dimensions. Then the N x M 
matrix 



built from gcd(M, N) copies of F as its diagonal and filled with zeros elsewhere, is the 
synthesis matrix of an M -element unit norm tight frame in C . Note that this is the sparsest 
unit norm tight frame one can construct by using altered DFT matrices as building blocks. 

It is this construction that we refer to as TDFTST in case that the number of vectors and 
the dimension are not relatively prime. 

3.2. General spectral tetris frames. We are now dropping any conditions on the sequence 
of prescribed eigenvalues, besides the frame condition, implying that the spectrum is positive, 
and the trace condition, implying that the sum of the eigenvalues has to equal the number 
of frame vectors of the unit norm frame to be constructed. We discuss a version of spectral 
tetris, which is capable of constructing a unit norm frame of M elements in with such 
spectrum, for any M > N. Note that now questions of how to order the eigenvalues come 
into play again and that therefore we will in general not get the sparsest frames possible 
from a spectral tetris like construction. The version of spectral tetris for this most general 
set up is the algorithm DFTST presented in table DFTST is a combination of the ideas 
developed so far and we will analyze it in the proof of theorem 13.121 

Definition 3.11. Given positive reals Xi,...,Xl we call a matrix a general Dl block for 
Ai, . . . , Xl, if it is derived from the DFT matrix Fl by multiplying the entries of the j-th row 

of Fl by J y- for j — 2, . . . , L — 1 and the entries of the first and last row of Fl by 



L 

respectively \p^, where < c\ < X\ and < cl < Xl, and if it has normalized columns. 
We call ci and cl the first, respectively last correction factor. 

Theorem 3.12. Let M > N be positive integers and (A n )^L 1 C (0, oo) in decreasing order 
such that Y2n=i = M. Then the algorithm DFTST constructs an M-element unit norm 



frame with eigenvalues (A„ 



n=l 



Proof. We construct the desired synthesis matrix successively one vector at a time by filling 
it with general Dl blocks, similar as seen in example 13.21 only now each time we check for 
the minimal size of a general Dl block that can be inserted. We do so by determining the 
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set S in line 4) of the algorithm. Indeed, for L £ S there does exist a general Dl block for 
Xj, . . . , Aj + £_i, since the condition on L is equivalent to 

1 _ :_£ < A i+ L ~ 1 
^ L ~ L ' 

and thus implies that the choice of Aj as the first correction factor results in a last correction 
factor, which is uniquely determined by the normalization condition of the columns, which 
is smaller or equal to Aj+l_i. 

Since we desire sparsity as a key property of the synthesis matrix, we choose the smallest 
possible such L, i.e. the minimum of S, in case S is not empty. If this minimum is 1 we just 
insert a 1, or in other words a D\ block, i.e. a standard unit vector. This is done in lines 
13) to 16) of the algorithm. Otherwise we insert a general D m [ n s block, see lines 18) to 25). 

The algorithm terminates after one performance of the lines 6) to 11). In lines 6) to 11), 
the first N — j + 1 rows of a general DM-k+i block are being inserted. This is done in case 
5 = 0, i.e. if we can not insert a general block for any L £ {1, . . . , N — j + 1}, or if 
after inserting a general D min s block, we would not be able to continue a spectral tetris like 
construction due to the fact that we would now have more rows than columns left to fill in 
the to be constructed synthesis matrix. □ 

4. Fusion frames with prescribed eigenvalues and prescribed dimensions 

Let M > N be positive integers and (X n )^ =1 C (0, oo) such that J2n=i = M. Given a 
sequence of dimensions, we ask the question of whether and how we can find a spectral tetris 
fusion frame for ~R N whose subspaces have those prescribed dimensions and whose fusion 
frame operator has the eigenvalues (\ n )n=v 

To get started, consider the following example of integer eigenvalues: 

(\ n )Li = (4,4,3,3,2,2). 

Given this sequence of eigenvalues, the spectral tetris frame in M 6 consists only of standard 
unit vectors: 

S = {ei, ei, ei, e 1} e 2 , e 2 , e 2 , e 2 , e 3 , e 3 , e 3 , e 4 , e 4 , e 4 , e 5 , e 5 , e 6 , e 6 }. 

The question we are asking above now takes the following form. We want to partition S into 
sets of pairwise orthonormal vectors, i.e. each set of the partition should not contain more 
than one copy of any standard unit vector. What sizes can these sets have? We start by 
considering the partition S = U n=1 Pn, where 

Pi = {ei,e 2 ,e 3 ,e 4 ,e 5 ,e 6 }, 

Pi = {ei,e 2 ,e 3 ,e 4 ,e 5 ,e 6 }, 

^3 = {ei,e 2 ,e 3 ,e 4 }, 

Pi = {ei,e 2 }. 

The sets of this partition have sizes 6, 6, 4 and 2. To get a different partition we cannot take 
any vector from Pj and put it into Pj if % > j, since this would destroy the ort honor mality 
of the sets. We can on the other hand take certain vectors out of Pj and put them into Pj 
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if i < j, without destroying the orthonormality of the sets. By doing so, we can for example 
easily find a partition into orthonormal sets of sizes 6,5,4 and 3. But, it is not possible to 
find a partition into orthonormal sets of sizes 6, 6, 5 and 1. The sequence 6, 6, 4, 2 majorizes 
the sequences of sizes of orthonormal sets which we can partition S into. Let us recall the 
notion of majorization. 

Definition 4.1. Given a = (a n )^ =l G M. N , denote by a} G M. the vector obtained by 
rearranging the coordinates of a in decreasing order. If (a n )^ =1 , (b n )^ =1 G M. N , we say (a n )^ =1 
majorizes (& n )n=i> denoted by (a n ) h (b n ), if ]T™ =1 > ]T™ =1 for all m = 1, . . . , N - 1 

and J2n=i a n = En=l 6 n- 

We will also use the notion of majorization between tuples of different length, by agreeing 
to add zero entries to the shorter tuple, in order to have tuples of the same length. 

We can use the idea of the above example to construct spectral tetris fusion frames in the 
general case of non integer eigenvalues. As above, we will determine a sequence of numbers 
depending on the given eigenvalues (A n )^ =1 and check whether or not this sequence majorizes 
the given sequence of dimensions. Again, the sequence we are going to determine will be the 
sequence of dimensions of a certain fusion frame for WL N having the eigenvalues (A n )^ =1 . We 
now introduce this fusion frame. 

Definition 4.2. Let M > N be positive integers and let (A n )^ =1 C (0, oo) have the property 
that Yln=i ^ n = The fusion frame constructed by the algorithm RFF presented in table® 
is called the reference fusion frame for the eigenvalues (A n )^ =1 . 

An example of RFF is included in the remarks following the proof of Proposition 14.41 
We will now use the reference fusion frame for (A n )^ =1 to tackle the question, whether 
or not a fusion frame for C , with a certain fusion frame operator and certain dimensions 
of the subspaces, is constructible via spectral tetris. In case it is constructible, the proof 
describes an algorithm to construct it. In the proof we will use the following notation. 

Definition 4.3. Let S be a set of vectors in <C N , and s 6 S. We say that a subset CCS 1 
is a chain in S starting at s, if s G S and the support of any element in S intersects the 
support of some other element of S. We say that C is a maximal chain in S starting at s if 
C is not a proper subset of any other chain in S starting at s. 

Theorem 4.4. Let M > N be positive integers, (A n )^ =1 C (0, oo) and let {dj)f =1 C N 
such that Y^n=i = Xln=i = M. Let (Ki)n=i be the reference fusion frame for (A n )^ =1 . 
If (dimV^) >: (d n ), then there exists a spectral tetris fusion frame (W n )^ =l for C N with 
dim W n = d n for n = 1, . . . , D and eigenvalues (A n ,)^ =1 . 

Proof. We show how to iteratively construct the desired fusion frame {W n )n=i m case the 
majorization condition holds. Let t and S±, . . . , S t be given by RFF for (\ n )n=v Let = Si 
for i = 1, . . . ,t. We add empty sets if necessary to obtain a collection (Wf)f =1 of D sets. 
If E*=i Wi \ - ck\ = then the sets (W?)fL x span the desired fusion frame. Otherwise, 
starting from {W^)f =l , we will construct the spanning sets of the desired fusion frame. 
Let 

m = max [j < D: dj ^ \Wj\} . 
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Note that YlT=i \ = Y^lLi^i by the choice of m, and YlT=i \Wf\ > YlT=i di by the 
majorization assumption. Therefore, d m > \W^\ and there exists 

k = max {j < m: \Wj\ > dj} . 

Notice that \W^\ < d m < d k < \W%\ implies \W^\ + 2 < |W fe °|. 

We now have to consider two cases. First, if there is some element w G W%, which has 
disjoint support from every element in W^, define {Wl)f =1 by 

f \ {w} iii — k, 
W l = pluW if i = m, (4) 
(W? else. 

Second, suppose there is no such element in W®. Partition U into maximal chains 
Ci, . . . , C r , say. Note that for i = 1, . . . , r, the sizes of the sets Cj PI and Ci fl differ 
by at most one, since, given v k G W® and t> m G W 7 ^, we know that v k and v m either have 
disjoint support, or their support sets have intersection of size one. Since \W^\ + 2 < \W%\, 
there is a maximal chain Cj that contains one element more from W® than from W^. Define 
{Wl)f =1 by 

((w»vjc 3 )\(c 3 nw») if i = A, 

^= WuC^^-nW:) ifi = m, (5) 
[w? else. 

In both of the above cases @ and (jSJ), we have defined (W/)^ such that 

D D 

Y^wi\ - d^ <y,w°\ - d t \. 
i=i i=i 

Note that {Wl)f =l satisfies the majorization condition in the sense that (IW 7 ^)) b (d n ). 
Thus if the sets of (W})? =1 do not span the desired fusion frame, we can repeat the above 
procedure with (Wl)? =1 instead of (W?)^ and get (Wf)f =l such that Ylf=i \\ W i\ ~ d i\ < 
J2f=i \ \ — di\. Continuing in this fashion we will, say after repeating the process / times, 
arrive at (Wj)f =1 such that J2i=i — di\ = 0, i.e. the sets of (Wl)f =1 span the desired 
fusion frame (W n )® =1 . □ 

Intuition suggests that for a given choice of eigenvalues, the dimensions derived from RFF 
for these eigenvalues in blockwise order, majorize the dimensions derived from RFF for the 
same eigenvalues in non-blockwise order. We do not investigate in this direction, since even 
two different blockwise orderings of the eigenvalues will in general lead to different sequences 
of dimensions of the reference fusion frame, as the following example shows. Given the 
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eigenvalues §, y , , STC produces the synthesis matrix 



1 1 
o o 
o o 







1 ^ ± 







^ 

12 

J_ 

12 



and thus the reference fusion frame for (§, y> if) * s 

V = (span{/i, / 5 , / 8 }, span{/ 2 , / 6 }, span{/ 3 }, span{/ 4 }, span{/ 7 }). 



Running RFF for a different blockwise ordering of the same eigenvalues, say y, ~, yields 
the synthesis matrix 



1 1 1 



o o J\ 
V 





•10 5 13 > 










1 



IS 



and thus the reference fusion frame for v 

W = (span{#!, g 6 }, span{# 2 , g 7 }, span{# 3 , g 8 }, span{# 4 }, span{# 5 }). 

Note that RFF is in both cases performed for a blockwise ordering of the given eigenvalues, 
yet the sequences of dimensions of the subspaces of the reference fusion frames V and W are 
different. 

In the case of constructing tight spectral tetris frames, questions of how to order the 
eigenvalues before performing the algorithm do not arise. Nevertheless, for the construction 
of tight spectral tetris frames of redundancy less than two, the majorization condition pre- 
sented in theorem 14.41 is not necessary either. Consider for example the case of a 10-element 
unit norm tight frame in C 7 . The synthesis matrix TDFTST constructs, is composed of 
two blocks D 2 followed by two blocks D 3 and the reference fusion frame derived from it has 
dimensions (2,2,2,2, 1,1). We might on the other hand construct a synthesis matrix for a 
10-element unit norm tight frame in C 7 by building its synthesis matrix using a D 2 block, 
followed by a D3 block, followed by another D 2 and D3 block. The reference fusion frame 
derived from this synthesis matrix has dimensions (2, 2, 2, 2, 2). 

We now show, that for constructing tight fusion frames of redundancy greater or equal 
to two, the majorization condition is also sufficient. In this situation, again, questions of 
how to order the eigenvalues do not arise. Moreover, the vectors in the spanning sets that 
RFF produces, are either standard unit vectors or linear combinations of two consecutive 
standard unit vectors. 

Definition 4.5. We call any standard unit vector G M. N a singleton. We call a linear 
combination aei + be^i of two consecutive singletons a doubleton and denote it by e[. 



By the structure of the synthesis matrix that spectral tetris produces and the fact that 
RFF sweeps through the synthesis matrix from left to right in order to fill the spanning 
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sets of the reference fusion frame, we see that RFF picks singletons first, i.e. if for some 
i — 1, . . . , N — 1 we have £ S' ril and e[ £ S n2 , then n\ < n 2 . 

Theorem 4.6. Let M > 2N be positive integers and (d n )^ =1 C N such that Y2n=i &n = M . 
Let (V n y n=1 be the reference fusion frame for (A„)^ =1 = (M, . . . , ^). Then there exists a 
tight spectral tetris fusion frame (W n )^ =l for 1R with dim W n = d n for n = 1, . . . , D if and 
only if (dimK) h (d n ). 

Proof. By Proposition 14.41 it remains to show that the majorization condition is necessary. To 
show this, let (S n Y n=1 be the spanning sets of the reference fusion frame, and suppose (F n )^ =1 
is some other partition of the frame vectors STC generates for (X n )^ =1 into sets of orthogonal 
vectors. It suffices to define a sequence of spanning sets of fusion frames (Si) n L 1 , j = 1, . . . , r 
for some r £ N, such that (SXii = (F n )%*, (S r n )ti = and (|^|) h (|5*|) for 

all j = 0, . . . , r — 1. We show how to construct such a sequence successively, i.e. we describe 
how to construct (S^" 1-1 ),^ 1 from (S 3 l ) n ^ =1 . 

So let (S®)nLi = (F n )n=i, fix j and suppose (S®)®li, • • • , (S^) n i, 1 have already been con- 
structed. If (S 3 l ) n i 1 = (S n Y n=1 , then j — r and we are done. Thus, we may assume that 
there exists 

n = min{n < Dj : S 3 n ^ S n }. 

Let m be the minimal integer in {1, ... , N}, for which S no and differ, i.e. for which one 
of the following holds: 

(i) No vector in S no is supported at m , but some vector in S{ is supported at m . 

(ii) e mo £ S no , but e mo £" S 3 nQ . 

(ill) e 'm £ S n , but e 'm £ S l - 

In each of these three cases we will describe how to construct (S 3 l +1 ) n ^ 1 from (5^)^^, such 
that (|5^ +1 |) y (\Sl\) and S 3 ^ 1 and S no either both contain no vector supported at m or 
the same vector supported at m . Note that by iterating this procedure, we will thus, after 
a finite number of steps, arrive at some (S^)^ which is identical to (S n Y n=1 . We now go 
through the three cases. 

Case (i): This case can actually not occur, i.e. it is not possible that S no contains no 
vector supported at mo but that S 3 nQ does. Indeed, S 3 no cannot contain a doubleton e^ reo _ 1 by 
the minimality of m , and it cannot contain a doubleton e' mo or the singleton e mo , since then 
RFF would have also picked just the same doubleton, or singleton respectively, for S no . 

Case (ii): Note that by the minimality of m , the reason for this case cannot be that 
contains a doubleton e' mo _ l . Thus, there are two subcases which could have created case (ii). 

(ii-a) S 3 no does not contain any vector supported at m . 
(ii-b) S 3 m contains a doubleton e' mo . 

Case (ii-a): There must be some k > n , such that e mo £ S 3 k and we define (S{ +1 )n=i by 
letting 

{S„ U{e mo } ifi = no, 
^\{e mo } ifi = fc, 
S 3 n else. 
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Case (ii-b): Again, identify the k > no, for which e mo G S k . Now, one of three things can 
happen. 

First, if S k contains the singleton e mo+ i, define (S , ^ +1 ) n ^ 1 by letting 

{( S L \ WmJ) u { e m , e mo+ i} if i = n , 
( s t \ {e mo . e mo+1 }) U {e^J if i = k, (6) 
S£ else. 

(To shorten our notation, in what follows we will use the following terminology for how we 
constructed (S^ 1 )^ 1 in ©. We say we constructed (5 , ^ +1 ) n l H L 1 from (S^)^ by switching 
e' mQ with {e mo ,e mo+ i} between and S 3 k .) 

Second, if S k contains no vector supported at m , define (Sl +1 ) n ^ 1 1 by switching e' mo with 
e mo between and S J k . 

Third, S k contains a doubleton e' mo+1 . This third subcase again has three subcases. The 
first being that contains e mo+ 2. In this subcase, we construct (S , ^ +1 ) n ^ 1 by switching 
{ e m ' e m +2} with {e mo ,e' mo+1 } between S^ Q and S k . The second being that contains no 
vector supported at m + 2. In this subcase, we construct (S'^ +1 ) n i 1 1 by switching e' mo with 
{ e m ) e m +i} between and S*^. The third being that S 1 ^ contains a doubleton e' mQ+2 . This 
third subcase will again have three sub-subcases, namely that S k contains e mo+ 3, no vector 
supported at m + 3 or a doubleton e' mo+3 . In the first two of those sub-subcases we can 
again define (S , ^ +1 ) n i 1 1 by switching vectors in the above fashion and the third sub-subcase 
will create three new further cases. We can continue our argument successively, considering 
three new subcases. Eventually, we must arrive at a point where we can define (Sl +l ) n ^ L 1 
by switching vectors, since we only deal with a finite number of vectors. 

Case (iii): In this case it is not possible that S 3 n contains e mo , since RFF picks singletons 
first and thus we would have e mo G S no instead of e' mo G S no . By the minimality of m , it is 
also not possible for S 3 no to contain some e' mQ _ x . Hence, S° n contains no vector supported at 
rriQ. There are three subcases which could occur under this circumstances. 

(iii-a) Si does not contain any vector supported at m + 1. 
(iii-b) e mo+ i G Sl Q . 
(iii-c) e' mo+1 G Sl . 

Case (iii-a): Identify k > n such that e' mo G S J k and define (S^ 1 )^ 1 by letting 

f^o^O tfi = n , 
S! +1 =<Si \{e' mo } if i = A f 
[Si else. 

Case (iii-b): Identify k > n such that e' mo G S J k and define (S , ^ +1 )^i 1 1 by switching e mo+ i 
with e' mo between Sl o and S k . 

Case (iii-c): This subcase has three new subcases in the same fashion as above, two 
of which lead to the definition on (Si +1 ) n ^ l 1 by switching vectors and the third of which 
again leads to three new subcases. So again we can continue our argument successively, 
considering three new subcases potentially over and over again, until eventually we arrive at 
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a point where we can define (S^ 1 )^ 1 by switching vectors, since we only deal with a finite 
number of vectors. □ 



In the trivial case of integer eigenvalues we can run STC and make the following observa- 
tion. 

Corollary 4.7. If (\ n )% =1 C N and {di)f =l C N such that £^ =1 A n = £n=i d n = M E N, 
where M > N, then a spectral tetris fusion frame with eigenvalues (A n )^ =1 and dimen- 
sions (d n )n =1 exists if and only if (a n ) >z {d n ), where a n = max{r: A r > n} for n = 
1, . . . , maxj =wv Aj. 

Proof. Note that in this case RFF produces the output t = max i=li jA r A, and that dim Vi = a« 
for i = 1, . . . , t. □ 
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DFTST: 


DFT Spectral Tetris 


Parameters: 


• 


Dimension IVgN. 


• 


Number of frame elements MgI 


• 


Eigenvalues Ai > ■ • • > A^ > such that J2^=i An = 


Algorith 


m: 


1) 


Set k = 1. 


2) 


For j = 1, . . . , iV do 


3) 




Repeat 


4) 




S = {L e {1, . . . , N - j + 1}: L < Yl^' X A^. 


/ 




If S = or if M - k = N - j and ^jf 115 " 1 A< £ N then 


} 




L = M - k + 1. 


7) 






8) 




For r = k, . . . , k -\- L — 1 do 


9) 






10) 




end. 


11) 




A; = A_j + i = • • • = Atv = U. 


12) 




else 


13) 




If min S = 1 then 


14) 




/* = ej. 


15) 




fc = Jfc + l. 


161 

/ 




A i = A j ~~ !• 


17) 




else 


18) 




L = min 5. 


19) 




(j = exp(^p). 


20) 




For r — k, . . . ,k + L — 1 do 


21) 










1 /i 1 ^+^2 \ , ,(L-l)(r-fc)„ 


22) 




end. 


23) 




A; = k + L. 


24) 




Aj = A J+1 = ■ ■ ■ = Xj + L_ 2 = 0. 


25) 




Aj+L-i = YJ l= j A; - L. 


26) 




end. 


27) 




until Xj = 0. 


28) 


end. 


Output: 




• 


Unit norm frame (f m )m=i C C . 



Table 5. The DFTST algorithm for constructing a unit norm frame with 
prescribed spectrum in (0, oo). 
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RFF: Reference Fusion Frame 








Parameters: 








• Eigenvalues (A n )^ =1 C (0, oo) such that Yln= 


i A n — 


M e 


N. 


Algorithm: 








1) Use appropriate spectral tetns algorithm lor 




=i to j 


let frame 


F = [Jm)m=l- 








2) t = maximal support size ol the rows ol b . 








o\ C. _ (7) f nr „• _ 1 f 








4) fc = 0. 








5) Repeat 








OJ fc = /C + 1. 








7) j = min{l < r < £ : supp fk D supp / s = 


for af 




S r }- 


8) $ = .s, u {./,), 








9) until k = M. 








Output: 








• Fusion frame (V^)* =1 , where = spanS'j for ? 


= 1,- 


..,t. 





Table 6. The RFF algorithm for constructing the reference fusion frame. 



